Search

The well known phenomenon of phase separation in synthetic polymers and proteins has become a major topic in biophysics because it has been invoked as a mechanism of compartment formation in cells, without the need for membranes. Most of the coacervates (or condensates) are composed of Intrinsically Disordered Proteins (IDPs) or regions that are structureless, often in interaction with RNA and DNA. One of the more intriguing IDPs is the 526-residue RNA binding protein, Fused In Sarcoma (FUS), whose monomer conformations and condensates exhibit unusual behavior that are sensitive to solution conditions. By focussing principally on the N-terminus low complexity domain (FUS-LC comprising residues 1-214) and other truncations, we rationalize the findings of solid state NMR experiments, which show that FUS-LC adopts a non-polymorphic fibril (core-1) involving residues 39-95, flanked by fuzzy coats on both the N- and C- terminal ends. An alternate structure (core-2), whose free energy is comparable to core-1, emerges only in the truncated construct (residues 110-214). Both core-1 and core-2 fibrils are stabilized by a Tyrosine ladder as well as hydrophilic interactions. The morphologies (gels, fibrils, and glass-like behavior) adopted by FUS seem to vary greatly, depending on the experimental conditions. The effect of phosphorylation is site specific and affects the stability of the fibril depending on the sites that are phosphorylated. Many of the peculiarities associated with FUS may also be shared by other IDPs, such as TDP43 and hnRNPA2. We outline a number of problems for which there is no clear molecular understanding.

The spontaneous formation of micelles in aqueous solutions is governed by the amphipathic nature of surfactants and is practically interesting due to the regular use of micelles as membrane mimics, for the characterization of protein structure, and for drug design and delivery. We performed a systematic characterization of the finite-size effect observed in single-component dodecylphosphocholine (DPC) micelles with the coarse-grained MARTINI model. Of multiple coarse-grained solvent models investigated using large system sizes, the non-polarizable solvent model was found to most-accurately reproduce SANS spectra of 100 mM DPC in aqueous solution. We systematically investigated the finite-size effect at constant 100 mM concentration in 23 systems of sizes 40 to 150 DPC, confirming the finite-size effect to manifest as an oscillation in the mean micelle aggregation number about the thermodynamic aggregation number as the system size increases, mostly diminishing once the system supports formation of three micelles. The accuracy of employing a multiscale simulation approach to avoid finite-size effects in the micelle size distribution and SANS spectra using MARTINI and CHARMM36 was explored using multiple long timescale 500-DPC coarse-grained simulations which were back-mapped to CHARMM36 all-atom systems. It was found that the MARTINI model generally occupies more volume than the all-atom model, leading to the formation of micelles that are of a reasonable radius of gyration, but are smaller in aggregation number. The systematic characterization of the finite-size effect and exploration of multiscale modeling presented in this work provides guidance for the accurate modeling of micelles in simulations., Comment: 19 pages, 5 figures, 2 tables, Corresponding Author: John E. Straub

Normal micelle aggregates of amphiphilic surfactant in aqueous solvent are formed by a process of entropically driven self-assembly. The self-assembly of reverse micelles from amphiphilic surfactant in non-polar solvent in the presence of water is considered to be an enthalpically driven process. While the formation of normal and reverse surfactant micelles has been well characterized in theory and experiment, the nature of dry micelle formation, from amphiphilic surfactant in non-polar solvent in the absence of water, is poorly understood. In this study, a theory of dry reverse micelle formation is developed. Variation in free energy during micelle assembly is derived for the specific case of AOT surfactant in isooctane solvent using atomistic molecular dynamics simulation analyzed using the energy representation method. The existence and thermodynamic stability of dry reverse micelles of limited size are confirmed. The abrupt occurence of monodisperse aggregates is a clear signature a critical micelle concentration, commonly observed in the formation of normal surfactant micelles. The morphology of large dry micelles provides insight into the nature of the thermodynamic driving forces stabilizing the formation of the surfactant aggregates. Overall, this study provides detailed insight into the structure and stability of dry reverse micelles assembly in non-polar solvent., Comment: 34 pages, 6 figure. In Ver2, GAP joined and re-plot figures, and we changed the title

Protein aggregation, linked to many of diseases, is initiated when monomers access rogue conformations that are poised to form amyloid fibrils. We show, using simulations of src SH3 domain, that mechanical force enhances the population of the aggregation prone ($N^*$) states, which are rarely populated under force free native conditions, but are encoded in the spectrum of native fluctuations. The folding phase diagrams of SH3 as a function of denaturant concentration ($[C]$), mechanical force ($f$), and temperature exhibit an apparent two-state behavior, without revealing the presence of the elusive $N^*$ states. Interestingly, the phase boundaries separating the folded and unfolded states at all [C] and $f$ fall on a master curve, which can can be quantitatively described using an analogy to superconductors in a magnetic field. The free energy profiles as a function of the molecular extension ($R$), which are accessible in pulling experiments, ($R$), reveal the presence of a native-like $N^*$ with a disordered solvent-exposed amino terminal $\beta$-strand. The structure of the $N^*$ state is identical to that found in Fyn SH3 by NMR dispersion experiments. We show that the time scale for fibril formation can be estimated from the population of the $N^*$ state, determined by the free energy gap separating the native structure and the $N^*$ state, a finding that can be used to assess fibril forming tendencies of proteins. The structures of the $N^*$ state are used to show that oligomer formation and likely route to fibrils occur by a domain-swap mechanism in SH3 domain., Comment: 12 pages, 8 figures, 9 supplementary figures (on 5 more pages), 2 supplementary movies (on youtube)

In this perspective we describe the critical role membranes play in modulating the structures of the Amyloid Precursor Proteins to produce the peptides involved in the Alzheimer's disease. Some of the key concepts related to protein aggregation including the potential role of the excited states of monomers in initiating protein aggregation are described., Comment: 1 Figure

To understand the mechanism of trimethylamine N-oxide (TMAO) induced stabilization of folded protein states, we systematically investigated the action of TMAO on several model dipeptides (Leucine, L2, Serine, S2, Glutamine, Q2, Lysine, K2, and Glycine, G2) in order to elucidate the effect of residue-specific TMAO interactions on small fragments of solvent-exposed conformations of the denatured states of proteins. We find that TMAO preferentially hydrogen bonds with the exposed dipeptide backbone, but generally not with nonpolar or polar side chains. However, interactions with the positively charged Lys are substantially greater than with the backbone. The dipeptide G2, is a useful model of pure amide backbone, interacts with TMAO by forming a hydrogen bond between the amide nitrogen and the oxygen in TMAO. In contrast, TMAO is depleted from the protein backbone in the hexapeptide G6, which shows that the length of the polypeptide chain is relevant in aqueous TMAO solutions. These simulations lead to the hypothesis that TMAO-induced stabilization of proteins and peptides is a consequence of depletion of the solute from the protein surface provided intramolecular interactions are more favorable than those between TMAO and the backbone. To test our hypothesis we performed additional simulations of the action of TMAO on an intrinsically disordered A{\beta}16-22 (KLVFFAE) monomer. In the absence of TMAO A{\beta}16-22 is a disordered random coil. However, in aqueous TMAO solution A{\beta}16-22 monomer samples compact conformations. A transition from random coil to {\alpha}-helical secondary structure is observed at high TMAO concentrations. Our work highlights the potential similarities between the action of TMAO on long polypeptide chains and entropic stabilization of proteins in a crowded environment due to excluded volume interactions. In this sense TMAO is a nano-crowding particle., Comment: 27 pages, 10 figures

A variety of neurodegenerative diseases are associated with the formation of amyloid plaques. Our incomplete understanding of this process underscores the need to decipher the principles governing protein aggregation. Most experimental and simulation studies have been interpreted largely from the perspective of proteins: the role of solvent has been relatively overlooked. In this Account, we provide a perspective on how interactions with water affect folding landscapes of A$\beta$ monomers, A$\beta_{16-22}$ oligomer formation, and protofilament formation in a Sup35 peptide. Simulations show that the formation of aggregation-prone structures (N$^*$) similar to the structure in the fibril requires overcoming high desolvation barrier. The mechanism of protofilament formation in a polar Sup35 peptide fragment illustrates that water dramatically slows down self-assembly. Release of water trapped in the pores as water wires creates protofilament with a dry interface. Similarly, one of the main driving force for addition of a solvated monomer to a preformed fibril is the entropy gain of released water. We conclude by postulating that two-step model for protein crystallization must also hold for higher order amyloid structure formation starting from N$^*$. Multiple N$^*$ structures with varying water content results in a number of distinct water-laden polymorphic structures. In predominantly hydrophobic sequences, water accelerates fibril formation. In contrast, water-stabilized metastable intermediates dramatically slow down fibril growth rates in hydrophilic sequences., Comment: 27 pages, 4 figures; Accounts of Chemical Research, 2011

Understanding the influence of macromolecular crowding and nanoparticles on the formation of in-register $\beta$-sheets, the primary structural component of amyloid fibrils, is a first step towards describing \emph{in vivo} protein aggregation and interactions between synthetic materials and proteins. Using all atom molecular simulations in implicit solvent we illustrate the effects of nanoparticle size, shape, and volume fraction on oligomer formation of an amyloidogenic peptide from the transthyretin protein. Surprisingly, we find that inert spherical crowding particles destabilize in-register $\beta$-sheets formed by dimers while stabilizing $\beta$-sheets comprised of trimers and tetramers. As the radius of the nanoparticle increases crowding effects decrease, implying smaller crowding particles have the largest influence on the earliest amyloid species. We explain these results using a theory based on the depletion effect. Finally, we show that spherocylindrical crowders destabilize the ordered $\beta$-sheet dimer to a greater extent than spherical crowders, which underscores the influence of nanoparticle shape on protein aggregation.

Energy transfer (relaxation) phenomena are ubiquitous in nature. At a macroscopic level, the phenomenological theory of heat (Fourier law) successfully describes heat transfer and energy flow. However, its microscopic origin is still under debate. This is because the phenomena can contain many-body, multi-scale, nonequilibrium, and even quantum mechanical aspects, which present significant challenges to theories addressing energy transfer phenomena in physics, chemistry and biology. In this paper, we describe our recent theoretical attempts to treat vibrational energy relaxation (VER) in biomolecular systems, including peptide, porphyrin, and protein. First we construct a reduced model using (instantaneous) normal mode analysis, and further add anharmonic coupling between vibrational modes. Using such a model combined with time-dependent perturbation theory for the anharmonic coupling, we derive non-Markovian VER formulas, which are applicable to small timescales such as sub picoseconds. We apply the VER formulas with full atomistic detail to various biomolecular systems, discuss the VER timescales and mechanisms, and compare with experiment if possible. We finally mention further theoretical prospects of the energy transfer phenomena in physics, chemistry, and biology., Comment: To be published in Adv. Chem. Phys. (2010)

Employing extensive quantum-chemical calculations at the DFT/B3LYP and MP2 level, a quartic force field of isolated N-methylacetamide is constructed. Taking into account 24 vibrational degrees of freedom, the model is employed to perform numerically exact vibrational configuration interaction calculations of the vibrational energy relaxation of the amide I mode. It is found that the energy transfer pathways may sensitively depend on details of the theoretical description. Moreover, the exact reference calculations were used to study the applicability and accuracy of (i) the quasiclassical trajectory method, (ii) time-dependent second-order perturbation theory, and (iii) the instantaneous normal mode description of frequency fluctuations. Based on the results, several strategies to describe vibrational energy relaxation in biomolecular systems are discussed., Comment: 18 pages, 6 figures, submitted to J. Phys. Chem. A

Using a time-dependent perturbation theory, vibrational energy relaxation (VER) of isotopically labeled amide I modes in cytochrome c solvated with water is investigated. Contributions to the VER are decomposed into two contributions from the protein and water. The VER pathways are visualized using radial and angular excitation functions for resonant normal modes. Key differences of VER among different amide I modes are demonstrated, leading to a detailed picture of the spatial anisotropy of the VER. The results support the experimental observation that amide I modes in proteins relax with sub picosecond timescales, while the relaxation mechanism turns out to be sensitive to the environment of the amide I mode., Comment: 10 pages, 5 figures, to be published in J. Phys. Chem. B

Vibrational energy transfer of the amide I mode of N-methylacetamide (NMA) is studied theoretically using the vibrational configuration interaction method. A quartic force field of NMA is constructed at the B3LYP/6-31G+(d) level of theory and its accuarcy is checked by comparing the resulting anharmonic frequencies with available theoretical and experimental values. Quantum dynamics calculations for the amide I mode excitation clarify the dominant energy transfer pathways, which sensitively depend on the anharmonic couplings among vibrational modes. A ratio of the anharmonic coupling to the frequency mismatch is employed to predict and interpret the dominant energy flow pathways., Comment: 13 pages, 4 figures, 2 tables, to be published in Chem. Phys. Lett

The need to understand the assembly kinetics of fibril formation has become urgent because of the realization that soluble oligomers of amyloidogenic peptides may be even more neurotoxic than the end product, namely, the amyloid fibrils. In order to fully understand the routes to fibril formation one has to characterize the major species in the assembly pathways. The characterization of the energetics and dynamics of oligomers (dimers, trimers etc) is difficult using experiments alone because they undergo large conformational fluctuations. In this context, carefully planned molecular dynamics simulation studies, computations using coarse-grained models, and bioinformatic analysis have given considerable insights into the early events in the route to fibril formation. Here, we describe progress along this direction using examples taken largely from our own work. In this chapter, we focus on aspects of protein aggregation using Abeta-peptides and prion proteins as examples., Comment: 27 pages, 7 figures, book chapter

Without invoking the Markov approximation, we derive formulas for vibrational energy relaxation (VER) and dephasing for an anharmonic system oscillator using a time-dependent perturbation theory. The system-bath Hamiltonian contains more than the third order coupling terms since we take a normal mode picture as a zeroth order approximation. When we invoke the Markov approximation, our theory reduces to the Maradudin-Fein formula which is used to describe VER properties of glass and proteins. When the system anharmonicity and the renormalization effect due to the environment vanishes, our formulas reduce to those derived by Mikami and Okazaki invoking the path-integral influence functional method [J. Chem. Phys. 121 (2004) 10052]. We apply our formulas to VER of the amide I mode of a small amino-acide like molecule, N-methylacetamide, in heavy water., Comment: 16 pages, 5 figures, 5 tables, submitted to J. Chem. Phys